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ABSTRACT 

We present CARMA 1.3 mm continuum data of the disk surrounding the 
young brown dwarf 2MASS J044427+2512 in the Taurus molecular cloud. The 
high angular resolution of the CARMA observations (0.16") allows us to spatially 
resolve for the first time the thermal emission from dust around a brown dwarf. 
We analyze the interferometric visibilities and constrain the disk outer radius 
adopting disk models with power-law radial profiles of the dust surface density. 
In the case of a power-law index < 1, we obtain a disk radius in the range of 
about 15-30 AU, while larger disks are inferred for steeper radial profiles. By 
combining this information on the disk spatial extent with the sub-mm spectral 
index of this source we find conclusive evidence for mm-sized grains, or larger, in 
this brown dwarf disk. We discuss the implications of our results on the models 
of dust evolution in proto-planetary disks and brown dwarf formation. 



Subject headings: circumstellar matter — brown dwarfs — stars: individual 
(2MASS J044427+2512) — planets and satellites: formation — submillimeter: 
stars 
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Introduction 



The investigation of disks around brown dwarfs (BDs) allows to test the physics of 
disk evolution and planet formation in physical environme nts which are very different from 



those characteristics of disks around young solar-like stars ( ILuhman et al.ll2007al : lApai et al. 



20081 ) . In the case of pre-main sequence (PMS) stars, disks are formed as a consequence of 
the conservation of angular momentum during the gravitational collapse of the protostellar 
envelope. Most of the mass forming a newborn star/BD comes from accretion through a 
circumstellar disk and characterizing the structure of accretion di sks around you ng BDs has 
the potential to constrain the mechanisms of their formation (see iLuhmanl 120121 ). 



Disk masses have been measured for a couple of tens of young BDs through far-IR 
and sub-millimeter photometry under t he assumption of optically thin emiss i on from dust 



and a ISM-like dust-to- gas mass ratio (IKlein et al. 



2003 



Scholz et all 120061 : lHarvey et al 



20121 : iRicci et al.ll2012bl ). However, because of the degeneracy between disk radius and disk 
surface density and geometry, photometry alone does not provide strong constraints on the 
radial distribution of mass in the disk. So far, the only direct obser vational estima t e of th e 
radius of a disk around a single young BD has been obtained by ILuhman et al.l ( j2007bl ). 
They combined Spitzer spectroscopy and high-angular resolution imaging with the Hubble 
Space Telescope for the nearly edge-on disk around 2MASS J04381+2611, and inferred a 
disk radius of about 20 — 40 AU from the optical images of the disk in scattered light. Also, 



Ricci et al.l ( 12012bl ) detected CO molecular gas from the disk surrounding the BD p— Oph 
102 with ALMA at 0.87 mm. They derived a range of possible disk radii: an upper limit of 
« 30 AU from the non resolved dust emission and a lower limit of « 15 AU from the CO 
brightness temperature, assuming that this emission is optically thick. 

In this Letter, we present the first spatially resolved observations of dust thermal 
emission from a young BD disk, 2MASS J04442 7+2512 (hereafter 2M0444). This object 
0.05 M , M7.25-spectral type, Luhmanl 2004), located in Taurus 1 Myr), is know n 
to harbor a disk from its IR excess knapp et alibooi iGuieu et al.ll20(m iBouv et al Iboosh . 
We targeted this BD disk, the brightest in the 1.3 mm survey by lScholz et al. using: 
high-angular resolution (~ 0.16") interferometric observations with the Combined Array for 
Research in Millimeter- wave Astronomy (CARMA) at a wavelength of about 1.3 mm. 

Furthermore, Spitzer spectroscopy of this disk provided evidence for /xm-sized dust 
grains, indicating that the very first steps of grain growth have taken place in this system 
( IBouv et al.ll2008l ) . We use our observations and others at sub-mm wavelengths to investigate 
the presence of larger, mm-sized grains in the outer regions of the 2M0444 disk and test the 
models of dust evolution in protoplanetary disks. 
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2. Observations and data reduction 

2M0444 was observed with CARMA array configurations A, C and D, providing base- 
lines lengths between about 11 and 1678 m. Observations in D-configuration were obtained 
on October 15 and 18 2011. Receivers were tuned at a local oscillator (LO) frequency of 
227.5 GHz, correspondent to a wavelength of about 1.318 mm. We setup the CARMA cor- 
relator to have a bandwidth of about 6.8 GHz for the continuum, and to observe at the same 
time the CO ( J = 2 — 1) molecular line, but no detection was found for this line. We used 
the same receivers and correlator setup for observations in C-configuration on January 15 
2012. Observations in A-configuration were carried out on December 10 2011. The LO was 
set at a slightly different frequency, i.e. 224.0 GHz, but considering the spectral index of 
2M0444 (see below) the flux of the source is expected to vary by < 4% at the two different 
LO frequencies. This is well below the estimated uncertainty in the absolute flux calibration 
which is about 10% for CARMA. 

The frequency-dependent bandpass was calibrated by observing 3C84 and 3C454.3, 
while the nearby quasar 0510+180 was used to correct for the time-dependent atmospheric 
and instrumental effects. During the observations in A configuration, the quality of these 
time-dependent corrections were checked by observing the nearby quasar 0431+206, located 
approximately 5 degrees north of 2M0444. This allowed us to quantify to < 0.05" the level 
of atmospheric seeing during these observations. This is much smaller than the angular 
resolution of our observations, i.e. 0.16", showing that the lower visibilities values at longer 
baseline lengths (Figured]) cannot be explained by atmospheric decorrelation. 

Uranus was used to calibrate the absolute flux scale for observations in D and C con- 
figurations, whereas 3C84 was used for data in A configuration, where Uranus is highly 
resolved. Since 3C84 is known to be a variable source in the millimeter, to evaluate its flux 
we interpolated between the fluxes measured on December 7 and 15 with the SMA at about 
1.3 mm. We obtained a flux density of about 8.3 Jy for 3C84 on the day of our observations. 
For all the four epochs of our observations we shifted the phase cent er to compensate f or the 



proper motion of 2M0444 (/i a = —1.7 mas/yr, fi$ = —26.5 mas/vr. iRoeser et al.ll2010l ). 



The MIRIAD software package was used for visibilities calibration and imaging. Left 
panel of Figure [1] shows the 2M0444 map after combining all the data from different array 
configurations. The integrated flux from the disk is about 5.2+0.3 mJy, where the uncertainty 
reflects the rms noise on the map. 
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3. Analysis of the observed interferometric visibilities 



Right panel of Figure [T] shows the real parts of the measured visibilities vs the projected 
baseline lengths of our CARMA observations. The decrease at longer baseline lengths indi- 
cates that the dust continuum emission of the disk is spatially resolved. This sets a lower 
limit to the angular extent of the disk, namely the angular resolution of the observations. 
At the Taurus distance of 140 pc, this corresponds to about 22 AU in diameter or 11 AU in 
radius. 

More precise information on the radial distribution of dust particles can be obtained by 
analyzing the measured visibilities. We use the same method adopted by llsella et al.l ( 120091 ) 
to fit interferometric visibilities of proto-planetary disks. They used a Markov Chain Monte 
Carlo (MCMC) to run large collections of disk models, they calculated for each disk model 
the predicted visibilities at the (u, v)-points sampled by their observations, and then derived 
estimates for the parameters describing the radial distribution of dust in the disk. 

To analyze the interferometric visibilities of the 2M0444 BD disk, we used two- layer 



(surface+midplane) models o f flared disks heated by the radiation of the central BD ( jChiang fc Goldreich 
19971 : iDullemond et al.l 120011 ). For the (sub-)stellar properties of 2M0444 we adopt ed an ef- 
fective temperature T e ff « 2838 K, a luminosity Lrd 0.028 L M as cons trained by iLuhman 
(12004 ) . and a distance of 140 pc. By considering the lBaraffe et al.l (120031 ) evolutionary mod- 
els of BDs and an age of 1 Myr, we de rived a mas s Mrd ~ 50 M Jup . This is within 10% 
from the value of 45 M Jup inferred by iBouy et al.l (120081 ) . We calculated the dust opaci- 
ties k v by considering porous spherical grains mad e of astronomical silicates, c arbonaceous 



materials and wate r ice (optical constants from IWeingartner fc Draind l200ll : IZubko et al. 



1996 ^ Warrenl 1984, res pectively) and used a s i mplifie d version of the abundances suggested 
bv IPollack et~aD jl994h . following Sicci et al.l feoiOal H). For the grain size distribution we 
adopted a power-law n{a) oc a~ q with q = 3.0 and truncated between a minimum size of 0.1 



/im and a maximum size which sets the value of f3 (k v oc z/^, e.g. iRicci et al.ll2010al ). 



We parametrized the surface density of dust in the disk as a truncated power-law in 
radius = Sioau(^/10AU) _p from an inner radius R in to an outer radius i? Q ut- For the 

disk inner radius we took R in = 0.05 AU, which is about the sublimation radius given by the 
BD properties, but all the parameters considered in our analysis do not significantly depend 
on the exact value of this parameter. Disk inclination i (angle between disk rotation axis 
and line-of-sight) and position angle P.A. (angle between north direction and disk major 
axis) have been left free to vary in the fitting process. However, disk inclination and position 
angle were not well constrained by our analysis. 



For our analysis we first fixed the value of the power-law index p of the radial profile of 
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the dust surface density. We adopted values o f p = 0, 0.5, 1, 1.5, similar to what found for 
disks around PMS stars (IGuilloteau et al.ll201ll ). We then launched several MCMC chains to 
check the convergence toward the absolute minimum of the x 2 -function in the (Eioau? Rout, h 
P. A.) parameter space of the disk models. This method allows to efficiently probe the region 
of parameter space around the y 2 -absolute m inimum, and to derive the confidence intervals 
for the fitting parameters (jlsella et al.ll2009h . 

Figure [2] shows the 2D contour plots of the x 2 -function on the (Eioau? Rout) parameter 
space for the different values of p. For a given p, there is a degeneracy between Sioau and 
i?out : for larger disks the disk surface density at a given radius has to decrease to conserve 
the total flux density which is proportional to the disk mass. There is also a clear correlation 
between p and the disk outer radius of the best-fit models: for disks with steeper radial 
surface densities, i.e. greater values of p ) the disk has to be larger to conserve the total flux 
density, which is dominated by the outermost disk regions. Our observations do not have 
enough signal-to-noise ratio in the faint outer regions of the disk, and the different p- values 
could not be discriminated by our data. This also means that the confidence intervals on 
R out will depend on the adopted value of p. These intervals are listed in Table [TJ In the case 
of p = or 0.5, the disk outer radius is between about 15 and 27 AU (at la). For larger 
values of p ) larger disks are consistent with the data, and in the case of p = 1.5 only lower 
limits for R out could be derived. The key point we want to highlight here is that for all the 
reasonable values of p we obtained a robust lower limit for R out > 10 AU. 



3.1. Grain growth in the outer disk 



Multi-wavelength obser vations of young disks in th e sub-mm have been used to probe 
mm-sized grains in the disk (IBeckwith fc Sargentlll99ll ). In the case of optically thin emis- 
sion in the Rayleigh- Jeans part of the spectrum, the flux density of the disk is given by 



Table 1: Confidence intervals for the disk outer radii R out as constrained by our analysis. 
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F v oc z/ 2+/? where f3 is the power-law index of the dust opacity k v oc zA Deviations from 
the Ray leigh- Jeans regime can be accounted for by radiative transfer modeling of the dust 
emission, while the contribution from optically thick emission can be quantified via resolved 
images of the disk ( ITesti et al.l 120031 ). Whereas tiny sub-micrometer sized grains are charac- 
terized by relatively steep sub-mm SEDs with j3 ~ 1.5 — 2, shallo wer spectra with 3^0 — 1 
have been typically found for young; disks around PMS star s (jAndrews fc Williams 2005 



Rodmann et al. 2006 : Lommen et al. 2007 : Ricci et al. 2010al Jbh . This result has been g< 



^en- 



(Draine 


2006; 


Natta et al. 


2007; 


Ricci et al. 


2012a) 



or larger, in the outer regions of the disk 



While grain growth to mm-sized particles have been found for several dozens of disks 
around PMS stars, the very faint sub-mm fluxes of BD disks make their characterization 
much more challenging. If compared to typical T Tauri disks, BD disks allow to probe the 
coagulation of solids to lower gas and dust densities which a re the most critical environments 
for the growth to mm-sized particles (IBirnstiel et al.ll2010l ). Also, models of dust dynamics 
predict that mm-sized particles should experience a faster inward radial drift in BD disks 
than in T Tauri disks (Pinilla et al., in prep.). Therefore, BD disks are ideal test beds for 
models of dust evolution in proto-planetary disks. 

So far, conclusive evidenc e for the presence of mm-sized grains in a BD disk has been 
obtained only for p Oph 102 (jRicci et al.ll2012bl ). They measured a sub-mm spectral index 
of 2.29 ± 0.16, and ruled out small (< 10 AU) optically thick disks from the analysis of 
the molecular CO emission. In the case of 2M0444, we measure a sub-mm spectral index of 
2.3 ±0.25 by combin ing our flux density measurement at 1.3 mm with 3.5 mm data obtained 
by iBouy et al.l ( 120081 ). i.e. 0.55 ± 0.12 mJy. Figure [3] shows the measured flux density and 
sub-mm spectral index overlaid with the predicted sub-mm fluxes for models of disks around 
a BD with the same properties as 2M0444. The sub-mm photometry of 2M0444 can be 
explained by two classes of models: relatively small optically thick disks with an outer disk 
radius R out ^10 AU, and larger, mostly optically thin disks with R out > 10 AU. In the 
former case, the disk emission is insensitive to both disk mass and dust opacity k v , and in 
fact all the possible values of f3 are consistent with th e data. This sho ws how photometry 
alone is not sufficient to constrain grain sizes in disks ( ITesti et al.l 120011 ). 



At the same time our CARMA observations spatially resolved the dust emission and 
ruled out models with R out ^10 AU: only the larger, mostly optically thin disk models 
with f3 = 0.50 ± 0.25 are consistent with both photometric and interferometric data for the 
2M0444 BD disk. This interval for (3 is consistent with most of the values constrained for 
proto-planetary disks around PMS stars, and provides evidence for the presence of particles 
with sizes > 1 mm in the outer regions of the 2M0444 BD disk. 
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Note that Figure [3] refers to disk models with a specific choice of the power law index 
of the dust surface density and disk inclination, i.e. p = 1 and i = 70 degrees, respectively. 
However, these parameters are not observationally constrained. We found that the results 
presented above are not affected by our ignorance on these parameters by investigating their 
effect in the ranges of p G [0, 1.5], i £ [0, 70] degree^]. 



3.2. Disk mass and outer radius: constraints on the mechanisms of BD 

formation 



The analysis outlined above shows that the dust emission of the 2M0444 disk is mostly 
optically thin. This means that a dust mass can be estimated from the sub-mm fluxes. The 
/3- values constrained by our analy sis, i.e. = 0.50=^ 0.25, correspond to dust opacities values 
of « I/} A=i.3mm ~ 0.2 — 7 cm 2 /gd us t ( IRicci et al.ll2010a[ ). We derive estimates for the dust mass 
of r>j 1 — 30 M Eart h, or ^ 0.3 — 10 M Jup in gas assuming an ISM-like gas-to-dust mass ratio 
of 100. These values correspond to a disk with a mass of ~ 0.6 — 20% of the mass of the 
BD, a nd are consistent with t he distribution of values for more massive disks around PMS 
stars (1 Williams fc Ciezall201ll ). 



The constrained values for the disk radius, i.e. R out w 15 — 30 AU for p — 0, 0.5 or R out > 
20 AU for p — 1, 1.5, are in line wit h previous estimates for BD disks detected with high 



angula r resolution optical ima ging ([Luhman et al.l l2007bl ). CO molecular line ( IRicci et al. 
2012b ). sub-mm photometry ( Scholz et al. 20061 ) . These values are very large for models 



of BD formation involving the ejection of pre-stellar objects in multiple systems, in which 
dynamical interaction during star-disk and disk-disk encounters signifi cantly modifies the 
appearance of BD disks ( e.g. Reipurth fc Clarke 2001 : Bate et al. 20031 ) . For example, the 
numerical simulations by Bate et al. ( j2003h show that young BDs formed in this way are 
typically surrounded by relatively small disks with radii w 10 AU, whereas only 5 — 10% of 
their BD disks should have R out > 20 AU. These numbers take in consideration the possible 
effect of viscous spreading which can lead to an increa se of the d i sk siz e on a timescale 
comparable with the disk age. Furthermore, as noted in iBate et al.l ( 120031 ). disks that were 
severely truncated and underwent following viscous spreading are expected to have a very 
low mass, because of both truncation and gas accretion during viscous evolution. This is 
not the case of 2M0444, whose disk-to-central object mass ratio inferred above is similar to 



1 The upper limit on the disk inclination comes from the low value of visual extinction of 2M0444, i.e. 
Ay « (|Luhmanll2004l ): a more inclined disk w ould absorb a fractio n of the sub-stellar radiation and would 
not be consistent with optical spectroscopy (see lSkemer et al.ll201lh . 
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what found for disks around PMS stars. 

Here we want to note that the sample of BD disks for which some information on the 
disk radius has been obtained so far is strongly biased toward the sub-mm bright objects. 
These might also correspond to the population of relatively large disks predicted by models 
of BD formation involving ejection processes, whereas the smaller disks with radii < 10 
AU would likely be fainter and therefore not detected at sub-mm wavelengths. A high 
sensitivity and angular resolution survey of homogeneously selected BD disks in the sub-mm, 
e.g. with ALMA, is neede d to constrain the predictions from these models on a statistical 



basis (jUmbreit et al.ll201ll ) and shed more light on the dominant channel for BD formation. 



4. Summary 

We present high angular resolution observations of the young BD disk 2M0444 obtained 
with CARMA at 1.3 mm. Our observations spatially resolve dust thermal emission for the 
first time for a disk around a BD. By analyzing the interferometric visibilities we infer disk 
radii of at least a few tens of AU for the 2M0444 disk. We find evidence for mm-sized 
grains, or larger, in the disk. Given the physical conditions in BD disks, i.e. lower densities 
and shorter timescales of solids drift than in typical T Tauri disks, this result provides the 
tightest test to models of dust evolution in disks, this topic being investigated in detail by 
Pinilla et al., in preparation. 
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Fig. 1. — Left) CARMA continuum map at 1.3 mm of the disk surrounding the young 
brown dwarf 2M0444. This map was obtained by combining all the data taken in A, C 
and D-configuration of CARMA. Contours are drawn at —3, 3, 6, 9, 12, 15a, where la = 
0.3 mJy/beam is the rms-noise. The yellow ellipse in the bottom left corner shows the 
synthesized beam, with FWHM = 0.36" x 0.44", PA = 80°, obtained with natural weighting. 
However, the longest projected baselines of our observations probe angular scales as small 
as 0.16". Right) Real part of measured visibilities for the 2M0444 disk vs projected baseline 
lengths. Points and errorbars represent the weighted means of the real parts and their errors, 
respectively, in the correspondent bins. 
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Fig. 2. — Contour plots of the x 2 -function projected into the (i? ou t, Sioau) plane, for different 
p-values. Contours were drawn at Ax 2 -values so that their projection on each axis gives the 
ID confidence intervals. These are drawn at the 68, 90, 95% confidence levels, going from 
the darkest to the lightest blue color, respectively. The best-fit models are represented by 
an orange star. 
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Fig. 3. — Sub-millimeter flux density vs spectral index for disk models for 2M0444. The 
black dot shows the sub-mm data for 2M0444. Each line represents the prediction of disk 
models with the same disk outer radius and dust opacity spectral index /3, as indicated on 
the plot, but increasing disk mass from left to right. Each model was computed by assuming 
a radial profile of the surface density E oc r _p , with p — 1, and a disk inclination of 70 
degrees. The adopted sub-stellar properties and dust opacities are discussed in Section [31 



